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O ' Abstract 

(N 

f^ I We address the problem of downlink beamformer design for signal-to-interference-plus-noise ratio 

•^ . (SINR) balancing in a multiuser multicell environment with imperfectly estimated channels at base 

'^ ■ stations (BSs). We first present a semidefinite program (SDP) based approximate solution to the prob- 

lem. Then, as our main contribution, by exploiting some properties of the robust counterpart of the 
optimization problem, we arrive at a second-order cone program (SOCP) based approximation of the 
^ , balancing problem. The advantages of the proposed SOCP-based design are twofold. First, it greatly 

reduces the computational complexity compared to the SDP-based method. Second, it applies to a wide 



range of uncertainty models. As a case study, we investigate the performance of proposed formulations 
when the base station is equipped with a massive antenna array. Numerical experiments are carried ou 
to confirm that the proposed robust designs achieve favorable results in scenarios of practical interest. 



> ■ 
QQ , when the base station is equipped with a massive antenna array. Numerical experiments are carried out 

en 

^T . Index Terms 

o 
en 

SINR balancing, massive MIMO, very large-scale antenna arrays, reduced complexity, interference 
channel, multicell beamforming. 



I. Introduction 

In practical wireless systems, it is virtually impossible to provide an error-free estimate of 
channel state information (CSI) to the transmitter. Although beamforming is very attractive from 
implementation and performance perspective, its applicability is reduced due to its sensitivity to 
channel estimation errors which may arise as a consequence of pilot contamination in multicell 
systems [HI, quantization effects due to digital processing [[2l etc. Motivated by this dilemma, 
various studies have been conducted to design 'uncertainty immune' precoders, see e.g., [|3l-[|6l 

The authors are with the Department of Communications Engineering and Centre for Wireless Communications, University 
of Oulu, Finland. Email: {mhanif, Itran, atoUi, markku.juntti)@ee. oulu.fi. 



and references therein. The key tool common to all the studies is the application of various 
important results from the robust optimization theory [0, [[8]|. For any optimization problem, 
the design of robust counterpart can potentially suffer from two major difficulties, namely, (i) 
hurdles in obtaining tractable representation of the robust counterpart of the original program 
thereby compelling to employ various approximations, and, (ii) once a tractable formulation is 
obtained an increase in the complexity of the robust counterpart is seen as compared to the 
original problem. This pattern is common to most robust designs pertaining to signal processing 
and communication applications in the literature. 

The significance of (ii) in designing uncertainty immune precoders is further enhanced when 
some of the parameters involved in the system setup take very large values. In this context, the 
recently envisaged large scale massive multiple-input multiple-output (MIMO) systems 0] can 
be considered. Indeed, such large-scale antenna arrays promise increased link reliability, better 
spectral efficiency and low power consumption at the cost of manyfold increase in the number of 
transmitter antennas compared with the traditional multiple-antenna systems. For instance, values 
of the order of hundreds of base station antennas have been proposed in flU. The sub-optimality of 
traditional precoding methods like zero-forcing, block-diagonalization is now well understood 
IfTOl . Any algorithm that is, for example, based on traditional mathematical programming is 
likely to outperform the heuristic approaches of zero-forcing etc. It is also pertinent to point 
out that the mathematical analysis of the present paper can be easily leveraged to the case of 
maximizing weighted sum rates based on the development presented in [fTOl . The traditional 
approaches of introducing robustness in the the precoder design can end up in a semidefinite 
program (SDP). The complexity of an SDP is highly sensitive to the precoder size (more details 
on this appear in Sec. IIII-BI) . and hence the SDP-based solutions can either incur appreciable 
computational cost or in certain circumstances the digital resources may not be sufficient to 
cater for the memory requirements of an SDP-based solution. On the other hand, second-order 
cone programs (SOCPs) are much more computationally efficient (again the details appear in 
Sec. IIII-BI) . and can certainly provide a viable alternative to designing algorithms for very large- 
scale antenna arrays. This motivates arriving at robust SOCP formulations for optimizing certain 
performance metric in modem communication systems. 

In this paper, we study the problem of signal-to-interference-plus-noise ratio (SINR) balancing 
in multicell multiple-input single-output (MISO) downlink or interfering broadcast channel with 



a realistic assumption of imperfect CSI. We focus on centralized base station (BS) control. 
The worst case design philosophy that is commonly employed in the existing literature is 
considered. We first show that the robust counterpart can be relaxed to an SDP, and, thus, 
can be (suboptimally) solved in conjunction with a bisection search. It appears that the SDP- 
based formulations provide a general solution to the robust design in many works, e.g., [JU- 
|l6l. However, the SDP-based methods may rely on a rank relaxation scheme, which is in 
general a suboptimal technique. Furthermore, the SDP-based approaches generally result in 
computationally expensive tractable robust counterparts. As our main contribution, we propose 
a robust design which is merely based on solving SOCPs, i.e., the proposed method does 
not represent much increase in complexity in comparison to the original version of the SINR 
balancing problem. This is accomplished by exploiting various properties of the constraints in the 
robust counterpart of the balancing problem. In particular, we avoid formulating the beamformer 
design by projecting it to the space of semidefinite matrices which normally results in a rank 
constrained SDP. More importantly, the proposed SOCP-based design can be used in a wide 
range of uncertainty models. We notice that the SDP-based design formulations commonly used 
in literature [|4l-[|6l are only applicable to the cases where the channel errors lie in an ellipsoid. 
As mentioned above, we also compare and contrast the SDP and SOCP solutions, particulary 
from the computational cost perspective, when the number of base station antennas is very large 
lUl. Finally, through numerical investigations, we show that the proposed SOCP-based solution 
offers comparable performance to the approach in [|6] and the SDP-based method when same 
uncertainty set (a ball) is used to represent channel perturbations. 

The rest of the paper is organized as follows. Section |ll] presents problem formulation, a 
solution for perfect CSI, and modeling of the balancing problem with imperfect CSI. Section HII] 
discusses in detail various solutions with imperfect CSI along with a comparative discussion about 
their properties. Finally, Sections |IV] and |V] describe numerical experiments and conclusions, 
respectively!!! 




Fig. 1. System model of a multicell MISO downlink channel. Green lines represent desired signals, while red ones denote 
interference. The serving base station for user k is denoted by bk- 



II. System Model 

Consider a system of B coordinated BSs and K users where each user is served by one BS. 
Each BS is equipped with T transmit antennas and each user with a single receive antenna. A 
sketch of the system model is presented in Fig. [B Interference originating outside the coordinated 
system is omitted. The serving BS for the A;th user is denoted by bk- The signal received by the 
kth user is 

K 

" (1) 



where h^^. ^ G C^^^ is the channel (row) vector from BS bk to user k, x/c G C^^^ is the 
transmitted signal vector from the BS b^ to user k and n^ ~ CM{0, a'^) represents circularly 



'We use bold lowercase letters to express vectors and bold uppercase letters to represent matrices. (.) , (.) and Tr(.) 
represent the Hermitian, transpose and the trace operators, respectively. C"'^ and R'"^'' represent the space of complex and 
real matrices (vectors) of dimensions given as superscripts, respectively. \M\ denotes the cardinality of set A4. [p]k represents 
the fcth component of vector p. \c\ and K(c) represent the absolute value and the real part of a complex number c, respectively. 
It denotes a T x T identity matrix. Finally, ||.||2 represents the I2 norm. 



symmetric zero mean complex Gaussian noise with variance o"^. The transmitted signal vector is 
defined as x^ = nikdk, where m^ G C^^^ is the unnormalized beamforming vector and dk is the 
normalized complex data symbol. The total power transmitted by BS b is Xlfcew ^^ (-^ [^fc^fe]) = 

II 1 1 2 

J2k(^u II^^*:||2' where the set Ub with size Kb = \V(b\ includes the indices of all users served by 
BS b. The SINR at user A;'s receiver is 

\K,k^k\^ .^. 

7fc = -^ (2) 

where the interference power in the denominator is divided into intra- and inter-cell interference 
components. 

A. Problem Statement and Solution for Perfect CSI 

For the case of perfect CSI, the maximin SINR balancing can be cast as 

maximize min a^Tfe G) 

mfc:EfeewJI™fclli<-Pi>.V6 k ^ 

where a^ are positive weighting factors. Using (|2]), we can equivalently reformulate (|3]) as 
maximize t 

subject to ||(hi,fcMi ■ ■ ■ hB,k^B (y^h < \AT^|hb„fcmfc|, V k, (4) 

||vec(Mfe)||2< Vn, V6. 
where M5 = [mz^^(i), . . . , uiu^(\Ui,\)\ includes the precoders of all users being served in the 6th 
cell and the operation vec(.) vectorizes the argument matrix by stacking columns. Furthermore, 
we can still find an optimal solution of (H]) even if h;,^ fciri/c, for all k, is forced to be real [|3l, 
[[TT|. [[T2II. In this way, the constraints in dU represent second order cone (SOC) constraints for 
fixed t. Therefore, the original problem can be solved as a series of SOCP feasibility problems 
using bisection search [[TTl . [fT3l . 

B. Modeling of Imperfect CSI 

In real systems it is impossible to achieve perfect transmitter CSI due to several reasons 
mentioned in, e.g., iH, JU. Hence, robust designs dealing with channel errors are of practical 
importance. We consider the channel uncertainty model in which the true channel vectors are 



of the form 



hfe,fe = hfe,fe + ^<5fe^fe[vi,,fc]i = hb,fc + Vb,fcAb,fc, V6, /c (5) 



i=l 



where hb^k represents the nominal (known) value of the channels, Ib^k ^ {I, 2, ... , T}, the vectors 
61 f^ (channel perturbation directions) form the rows of A^ ^ G C^^.k/xT ^^^ v^ fcAfe^ gives the error 
vector in the downlink channel from BS 6 to user /c Q, lfT4ll . We denote by S the uncertainty 
set that includes all channel error row vectors v^^. As seen in ^ the above model assumes 
that the uncertainty vector affects the data in an affine manner. This philosophy has been widely 
used, e.g., in [|5l, O etc. In addition to affecting the true channels in an affine manner, the error 
vectors are also constrained to lie in an uncertainty set St^k as 

<Sb,k = {vfe,fc : \\vb,k\\ < Pb,k,'^b, k} (6) 

where || ■ || is an appropriate norm specified by the parameter p^^, and is chosen based on how 
one wishes to model channel uncertainties. Normally, modeling channel errors by exploiting 
their statistical nature is prone to numerous difficulties. To name a few, it requires information 
about the statistics of the error vectors which is mostly not available on account of myriad of 
phenomenon involved in the channel estimation process. Then even if some information about 
the statistics of the error vectors is exploitable, a part from simple linear constraints contaminated 
with Gaussian errors, it is virtually impossible to arrive at exact tractable versions of stochastic 
constraints. Motivated by the dilemma, the uncertainty in channels is modeled by norm-bounded 
sets ^. With such modeling, it does not remain necessary to know information about the, say, 
probability law that the uncertainty vectors follow. Further to this, as we will see in the discussion 
to follow, norm-bounded uncertainty sets model various real world scenarios very well. One more 
advantage of such modeling is that in several cases of interest, the norm-bounded uncertainty 
models permit either exact tractable formulations or good approximations [|T4l. [[T5l . 



C. Worst Case Design Formulation 

We will concentrate on the worst case robust optimization approach of JH], [[T4l that has 
been traditionally used in the existing literature for different problems [|4l-[|6l. The worst case 
approach amounts to satisfying the constraints for all possible channel vectors. Hence, the robust 



counterpart of © is written as 
max. t 

mfc,t 

B 

s. t. —\hb^^kmk\ > ^ |hf,^,fcmi| + ^ ^\hb,kmi\ +a'^, W k, W {vb,kM,k ■ ^b,k ^ Sb,k} 

5^||mfc||2<P,,V6. 

(7) 
We note that the formulation in (|7]l is intrinsically intractable owing to its semi-infinite nature 
i.e., finite optimization variables and infinite constraints. 

Remark 1: It is worth mentioning here that for the case of receivers equipped with multiple 
antennas (and hence the possibility of transmitting multiple data streams), an option could be 
to employ receiver combining matrix and study the balancing problem on per stream basis. The 
problem is then exactly similar to our case for a given receiver processing matrix per user. Here 
we stress that even in the presence of perfect CSI, joint transmitter-receiver processing matrix 
design is a difficult nonconvex problem. 

III. Robust Beamformer Designs for Imperfect CSI 

A. SDP-Based Robust Design 

The solution to problem (|7]) significantly depends on the type of the uncertainty set. In the 
context of SDP design, the most commonly considered uncertainty model is the one where error 
vectors V5 ^ are bounded in a ball of radius p^ fc, i.e., Sb^k is expressed as 

Sb,k = {vfe,fc : ||vfe,fc||2 < pb,k,'^b,k}. (8) 

It is clear from (|5]) and ^ that hb^,. is assumed to lie in an ellipsoid centered at h^ fc, which 
is characterized by pb^k and A^fc. This type of channel error model is commonly known as 
ellipsoidal uncertainty model in the related literature. For practical channel estimation schemes 
it is known that the channel estimation error follows Gaussian distribution [|T6l . Most of the 
probability content of multi-dimensional Gaussian density is localized in its certain region. This 
clearly motivates modeling the error using an ellipsoid. Further, when vector quantization is 
used at the receiver, quantization errors can also be approximated by ellipsoids [fTTll . For this 
specific uncertainty set, we will show that the robust counterpart in © can be solved using SDP 
approximations. We note that in © the uncertain part of h^j. /,. varies in the same set on both 



sides of the inequality. It is well known that this renders the problem intractable 01, [lU, [fTSl . 
However, we will see that, after a suitable relaxation, this constraint can be written in a tractable 
form. To this end, we define P^ = mfcinf, Q5 = Y^keUt^k and W^ = ^P^ - Y.ieUt^\kPi- 
Then by using slack variables we observe that the constraints involving perturbed channels can be 
cast into tractable linear matrix inequalities (LMIs) using the so called 5-Lemma or 5-Procedure 
[HI, [fTSl . After some manipulations, (|7]) can be equivalently rewritten as 

maximize t 

bk,k h^ ^ ^ b„L b„k ^ 0, V A: (9a) 

K.fcWfeA" , V,fcWfch« ,, - Tk - Xb„kpl,kJ 

-Ab,kQbA.fk + Ab,fcl -Afe,fcQfeh^fc \ 

I ^0VA;,6 7^6fc (9b) 
—^b,kQbAbi, —hb^kQb\^k + ^6,fc ~ \kPb,kJ 

B 

J2 trace(Pfc) < Pb, V6, ^ tb,^ + f^^ < r^, r^ > V fc (9c) 

keUi b=l,b^bk 

Xb,k > V 6, fc, tb,k >0\/ k,b^ bk, Pfc ^ 0, rank(Pfc) = 1 V fc. (9d) 

Therefore, after ignoring the nonconvex rank constraints, bisection search over t can be used to 
obtain covariance matrices. However, we cannot guarantee optimality (see |fT9ll ) of the proposed 
solution if we obtain the rank of precoding matrices greater than unity. Similar rank relaxation 
approach was also adopted in the recent works, e.g., [51. We may need some randomization 
procedure [[201 to extract the beamformer nifc if rank(Pfc) > 1. Nonetheless, randomization trick 
may not always be useful [[5][. 

B. SOCP-Based Robust Scheme 

The robust counterpart of any optimization problem can potentially pose two issues, one 
related to its tractability, and the other related to its complexity. Very often, the worst case 
principle leads to an intractable problem, since, as noted previously, the robust counterpart is 
an optimization problem over an infinite set of constraints. Furthermore, commonly employed 
approximation schemes usually increase the complexity of the original problem by one degree, 
i.e., a linear program becomes an SOCP and an SOCP transforms to an SDR In what follows, 
we propose a robust design which is merely based on iteratively solving SOCPs, i.e., we attempt 



to minimize the complexity of the robust version of the balancing problem. Interestingly enough, 
we note that the SOCP-based scheme can also encompass a wide variety of uncertainty sets. In 
order to emphasize the capability of the SOCP scheme to handle a variety of uncertainty sets, 
we will not specify any particular norm to represent the uncertainty set. 

We will arrive at a reduced complexity tractable robust scheme by incorporating uncertainty 
and exploiting the structure of the SOC constraints in dU. To start with, we consider a relaxation 
of ([7]), which is written as 



maximize 



t 




hf.,k 



^h^,k + ^sl,,k[^h^,k] 



i=l 



nifc 



K^i.fc ••• ZB^k o-)^|L > 0,VA;, 



Vvfc, i. 



\^hk,k\\ ^ Pbk,k 



Zb.k' 



M^ 



l-b.k 



^K,k[^b,k]i 



(10a) 



^b,k 

i=l 

Vvfe^fc : \\^h,k\\ ^ Pb,k, 
Zb,k>0,Wb,k, ||vec(M,)||2< 
, nik G C'^''\zb,k G M,M, 



> 0,V6,fc, 



b, 



V6, 



(10b) 
(10c) 



where Ck = ^1 + ^, rrife G 0^'''^,Zb,k G M, M5 = [mu,{i) , . . . , m.u^QUb\)] are optimization 
variables and we have removed the absolute value, and only consider the real part of the left side 
of the constraints in (|10al) . Unlike the non-robust version of the problem dU, we cannot force the 
imaginary part of [h;,^. ^ + J2i=i ^l k[^bk,k]i]i^k to zero for all channel error realizations. Since 
for a complex number x, \x\ > ^(x), a feasible point for (llOal ) is also feasible for the exact robust 
counterpart given in (|7]). That is to say, to arrive at a reduced complexity approach, we consider 
a conservative approximation of the exact robust counterpart of the problem. For notational 
simplicity, when clear from context, we avoid mentioning the real operator ^(■) explicitly from 
this point onwards. 

Now we make a key manipulation by substituting v^ ^ = 96^ — <pf,k ^^ all b, k, in (|10bl) and 
3fJ(vfe„fc) = ^{eb„k) - ^{<Pb„k) such that ^{db^k) > and ^{ipk„k) > in m^. After this 



we obtain a relaxation of (flOl) 
maximize t 

subject to 



5R([c|. 



6fc,fcJ« 




- \\[Ziu 



ZB,k (T) 



ve,„,, cij,^,, : P(e,„fc) + 3?(tl56„.)ll < pL. 



^fe,fc- 



M« 






l2>0,Vfc, 
(11a) 



Kk + Y.^U^^^bA^-^'^ 



bMi) 



i=l 



> 0,V6,A;, 



V0b,fc, f^hk ■ \\\^b,k\ + \<Pb,k\\\ ^ Pb,k^ 



Zb,k>0,^b,k, ||vec(M,)||2< 



b, 



V6, 



(lib) 
(lie) 



where [Qb^.tli, [4>fefe,fe]i denote the zth components of Qbk,k, ^bk,k for all b, k, respectively, and for 
a vector y the symbol |y| represents that [|y|]j = |[y]j| for all i. Another change introduced in 
(fTTI) is that for all b, k we have replaced pb^^^k with p[,^. The motivation for this variation of the 
uncertainty set parameter p^^.t becomes clear as we outline the fact that splitting the uncertainty 
vector V{,j. fc into a difference of two vectors and manipulating the left side of uncertainty sets 
as done in (fTTj) transforms the problem into nearly a safe approximation of its original version 
while also rendering it tractability . 

Remark 2: As noted earlier, it appears difficult, if not impossible, to cast the worst case robust 
counterpart (|7]) into its exact equivalent tractable formulation. For example, in the first approach 
based on SDP formulation, we had to drop the unit rank constraints to arrive at a tractable 
representation. Naturally, to arrive at an SOCP representation of the problem, we have to resort 
to an approximation of the original feasible set i.e., 

Oorig = {Optimization variables in (|7]) such that all constraints in (|7]) are satisfied} (12) 

with its tractable subset that may also include some additional analysis variables. 
In our case, it should be a feasible set for an SOC problem. By doing so we can ensure that 
a solution for the approximation is definitely feasible for the original optimization program, 
and thus promises safety in a sense that we do not violate the original constraints. By just 
considering the real part of the left side of constraints in (|10al) . we follow this strategy. It is 



worthwhile to note that the approximation used in (111 al l may compromise safety. However, our 
numerical investigations in Sec. |IV] reveal that the safety of the proposed SOCP procedure is 
almost guaranteed. It is evident that by deriving a subset of the original feasible set, we may 
have an overly conservative approximation. Hence, manipulating the parameter p^ ^ as p^ ^ may 
provide some flexibility to overcome the conservativeness of the proposed approximation ^2\\ . 
Next, let us first focus on the set of constraints in (lllbl ) and rewrite it in a form similar to 
the one presented in (llObl l 



Zhk — 



M« 



Y^^lk 



Vfe,fc 



«=1 



> 0,V6, A;,Vv6,fe : ||v6,fc|| < p'b, 



(13) 



It is worthy making an important observation now. A set of optimization variables satisfies the 
constraints in (lllbl ) if and only if it satisfies the set of constraints in (IT3T ). For a given b, k, let us 



assume that Zb,k and M5 are infeasible in (lllbl) . i.e., there exist G^fc, cf)^^ and |||9bfc| 
p[ ^ such that 



Zb,k 



M" 



'6,fc 



V + $^<5t^fe([0^fc] 



[* 



b±H> 



i=l 



<0. 



\'^bA\\< 



(14) 



Let [vfe,fc]i = [0fe,fc], - [^b,k\i for aU i. Thus it is easy to see that \[^b,k\i\ < \[^b,kU + l[fl^fe,fc]il- 
Therefore, we obtain \\^b,k\\ < |||0fe,fc| + I4'6,fclll ^ P'bk^ ^'^d hence (fT3l) is also infeasible. Next, 
we assume conversely that for a given 6, k, Zb^k and M^ are infeasible in (ITJj ). i.e., 

'b,fc -1 H 



Zb,k- 



M? 



ife,fc 



i=l 



<0 



(15) 



for certain ^/b,k such that ||v6,fc|| < p^^,. Let [9^,^]* = (1 - '^b,k)[^bAi ^nd [^b,k]i = -'^b,fc[vfe,fe]i, 
where ?9b,fc G [0,1]- With this substitution, it is seen that \wb,k]i = [Qb,k]i — [<Pbk]i- Similarly, 
these substitutions imply \[Qb,k]i\ + \[^b,k]i\ = 1(1 - i^b,k)\\[vb,k]i\ + | - i^b,k\\[^b,k]i\ = |[vfe,fc]i|, 
and, thus, |||96_fe| + Itt^b./dll = ll'^b,fcll < p'bk- Therefore, the variables Zb^k and M;, are infeasible 
in (lllbl) as well. Hence, we conclude that the feasibility of the constraints in (|llbl) implies the 
feasibility of (fT3l) . and vice versa. 

The main goal of the development so far is to approximate (lllal) and (lllbl) by SOC constraints 
so that the resulting robust counterparts in (lllal i- dllcl i can be cast as an SOCP for fixed t. Since 
(lllal) and (lllbl) have the same form, it is sufficient to concentrate on tackling the more difficult 
set of constraints in (|llbl ). We use the concavity of the negative norm to bound (|llbl ) from 



below as 



Zhk- 



lM"b" II - 

|-'-^-'-b^6,fc||2 



Mi 



h.k 



Y.^U[Qb,kh-['p 



b,k\i) 



I- 4=1 



>o, 



yb,k,'ieb,k,<Pb,k- lll0fc,fcl + |tl^6,fclll < p'b,k- (16) 

Again using the concavity argument, the left side of constraints in (fT6l) can be further lower 



nded 


as 




Zb,k 




4=1 



IKkMbne 



bM 



\{si,kMbn-[(i>b, 



kit) 



> 0, V6, k 



^Qb,k,(i>b,k ■ \\\Qb,k\ + |tt»fe,fc||| < p'b,k- (17) 

Reading the inequalities from (fTTl) backwards, and recalling the equivalence of (|llb|) and (fT3l) . 
we observe that a solution of (flTl) is also feasible for (flST ) or (llObI) . 

Remark 3: Before presenting a tractable formulation of the constraints in (flTI) . we again note 
that the optimal solution of the proposed SOCP relaxation ideally should also be feasible for the 
original worst case robust counterpart in (0. Therefore, if Osocp represents the feasible set for 
the SOCP relaxation, then Osocp ^ Oorig should hold. This will imply both safety and tractability 
for the proposed SOCP approximation. Although the transformations that will lead to an SOCP 
formulation for (flTl l ensure both these factors, the same cannot be observed for the constraint in 
( lllal ). Nevertheless, as also noted above, we will see in Sec. |IV]that the relation Osocp ^ Oorig 
almost remains valid at least for the cases considered. Being nearly a subset of the original 
problem, the proposed approximation can be rather conservative, as also noted in ifTSl . [|2T1l . To 
provide more flexibility in this regard, and as mentioned above, we make pb^k a design parameter 
and replace it by p^ ^ in (fTTl) . With the introduction of this maneuver, we may be able to improve 
the achieved objective, albeit this may come at the cost of degradation in achieving it for given 
realizations of channel errors as we probe in the results section. 

Let us define 

-||{5^M,}«||^. (18) 



fi{Mb, Zb,k,iib,k 



Zb,k - \\Mfhl,\\^, f2{Mb,Sl,) 



We note that /2(Mb, dlf.) = /2(Mb, -(5^^). With the above definitions and using the fact that 



||/i;m||2 = |/i;|||m||2, the constraints in (fTTj ) can be equivalently written as 

ll|ei,.fc|+|*6,fclll<p;,,fc^ I ' J 

(19) 
The constraint in ([19] ) can be cast into tractable form using [ |l5i, Theorem 1], which is stated as: 



Theorem 1: Working in the real domain, given a function /(x, U) that is concave in data U 
for all given x and scales linearly with the data, we consider a constraint of the following form 

min /(x, U") + V [/(x, 5^Mui], + /(x, -¥)\M,\ > (20) 

ui,U2>0:||ui+U2]|<tJ ^ — ' 

3 

where U" is the nominal part of the data, S^ is a vector representing perturbation direction in 
the jth component of the data and ui and U2 are real vectors of appropriate dimensions and the 
norm in (|20l) satisfies the property [jTSl Eq. 6] 

l|u|| = |||u||| (21) 

where |u| = (|ui|, . . . , \ud\). The constraint d20l) admits an equivalent representation of the form 

/(x, U") > a;||y||*, where [y]j = max{— /(x, 5-'), — /(x, — 5-')} > and ||y||* = max s'^y is 

||s||<l 

the dual norm of y. 

Proof: The proof of the theorem is available in [[T5l Theorem 1]. However, for the sake of 
completeness and for demonstrating its applicability on (1191) it is relegated to the Appendix. ■ 

It should be emphasized that the norm in (|20|) can be arbitrary, as long as it satisfies (|2T]) . 
meaning that the proposed SOCP-based scheme presented next is applicable to a wide variety 
and combinations of norms and thus uncertainty sets. For the special case of I2 norm, the norms 
remain I2 because of the self dual property of the I2 norm. Following similar steps used to tackle 
the constraints in (|llbl ). we can easily see that the uncertain constraints in (|llal) can be cast in 
a form that is amenable to applying Theorem [1] However, some important observations should 
be re-stressed at this point. Although the constraints in (II lal) are linear, the solution of this 
approximation does not necessarily imply (llOal l. This differs from the previous scenario where 
the conversion of the constraints in (II Obi ) to (lllbl ) is safe. Therefore, obtaining a safe, tractable 
and least possible conservative version of (llOal ) is left as an open question for future research. 

With the aid of Theorem \T\ the approximate robust counterpart of the original problem can 



be written in the following tractable form 
maximize t 

T 



subject to Ckhb^,kmk - 



zik ■ ■ ■ ZBk o- 



> pL,kLb,,k,yk {lid.) 



Ckdl^^^mk + K,,fc]g > 0, V/c, g = 1, . . . , 4^,^ (22b) 

- CkS\^,,mk + [(\hkAq > 0,VA:,g = 1, . . • ,4,.,^, ||qb,,fc|| < ^fefc,fc,V/i; (22c) 

- ||{5^Mf,}«||2 + [^i,,,], > 0,z = 1 . ..h,k,yb,k, \\ii,J < Uh,k,yb,k (22e) 

||vec(Mfe)||2 < ^/Pt, Wb (22f) 

where nik G C^''\ M^ = [mu^{i), ..., mu,(\u,\)], Zb^k G K, L^^^fc G R, Ub,k G M, qb^,fc G C^"!^-", ^b,k ^ 
Qk,k are optimization variables. The above optimization problem represents a tractable approx- 
imation, in the form of SOCP in conjunction with bisection search, of the robust counterpart 
of the problem under consideration. In the following, we provide some remarks regarding the 
tractability and reduced complexity of the proposed SOCP-based robust design. 

Tractability: We emphasize that while the SDP-based solution is only applicable to ellipsoidal 
uncertainty models, the SOCP-based approach is flexible enough to deal with other types of 
uncertainty sets. For example, in certain situations, the errors in each of the individual terms of 
the channel vector are bounded i.e., |[vfofe]j| < ^bM for all b,k,i. This amounts to saying that 
||vfe,fc||oo < ^b,k- In fact, in practical systems where each entry of h;,;,, is quantized independently 
at the receiver and fed back to the corresponding transmitters, the interval uncertainty model 
is more appropriate [22]. Clearly, this uncertainty model can be easily handled with the above 
approach since the dual of the l^o norm is well known [|8l. In other situations, it may happen that 
the entries of the uncertainty vector are symmetrically random and bounded. In such scenarios it is 
well known that the perturbation set can be represented as the intersection of the I2 and /qo norms 
of V;, fc ll23l . For this uncertainty model, it may be difficult, if not impossible, to straightforwardly 
use worst case design philosophy. Hence, the SDP-based method is not applicable and problem 
^ appears to be intractable. However, the SOCP-based approach admits tractability in the 
approximate solution of the robust counterpart using the dual of the I2 fl /qo norm [|23]| . 

Complexity Reduction: The SOCP-based robust design also offers a great reduction in compu- 



tational complexity compared to the SDP-based method. In what follows, we give a complexity 
comparison of the SDP- and SOCP-based solutions for the special case where A5 ^ = I^ for all 
b, k, which is commonly considered in the related workso First, let us focus on the equivalent 
representation obtained using Theorem [1] and explore it by considering any robust equivalent 
constraint (without loss of generality) from (I22el ). We note that under this setting each entry of 
a channel can be written as [hb,k\m = [h-b,k]m + [Sb,klm[vb,k]'m, 1 < m < T, where v^^fe belongs 
to the uncertainty set defined in ([8]). The vector y^k corresponding to the equivalent formulation 
of Theorem [T] becomes 

ylk = [\[T^U,{l)]m[Sb,k]m\, ■ ■ ■ ,\[mu,(\Ut\)]m[Sb,k]m.\],^b,k. (23) 

With this type of y^ ^ it has been shown in ifTSl that for ellipsoidal uncertainty set, instead of 
having multiple additional constraints of the type mentioned above, we can stack all correspond- 
ing variables into one SOC constraint, || M^f, fc||2 < t^b,k, and one variable i'b,k for all b, k. Similarly, 
the constraints involving user {bk, k), shown in (I22bl) - (l22cl) . can be greatly simplified. 

To provide a complexity comparison, we base our discussion on the simplification noted above 
and focus on an arbitrary bisection step. According to (|22al i- (l22fl) . the number of real optimization 
variables per bisection iteration of the SOCP-based robust design is ATBK + 2BK + 2KT + K. 
More specifically, there are BK constraints of real dimension {2TKb + 1) that occur thrice 
including the power constraint. Again using the above mentioned simplification, we obtain two 
constraints of real dimensions B + 2 and T + 1 that are K in number. Combining all these the 
worst case per iteration complexity of the SOCP approach approximates as 0[Kb{TBK)^) [|24l|. 
ESU . The per iteration complexity of the SDP-based method is found to be 0{{KBTf) El, 
ll26l . which is clearly higher than the SOCP counterpart. Further to this, based on ll24l . [|25l . the 
worst case estimate of the number of iterations needed to arrive at a numerically acceptable value 



of the SOCP-based design is 0{\/KB). As similar calculation reveals that such an estimate for 
the SDP-based method results in a higher value of 0{-\/KTB) on account of its dependence on 
the size of the matrix inequalities. A more detailed exploration that compares run times of the 
proposed approaches with different solvers is given in Sec. ITVl 

The same arguments in this part also apply to case where the entries of channel vectors undergo independent perturbations. 



IV. Numerical Results 

In order to compare the performance of the proposed approaches we report results of numerical 
simulations in this section. For all simulation setups, we consider a system of two cells (B = 2), 
while the number of users per cell is mentioned for individual numerical experiments. The 
channel vector from the BS b to user k is given by h;, fc = y^b^hfe ^ where Kb,k represents both 
the path loss and the shadow fading and h^^ follows CJ\f{0, 1). In Figs. [2l|4]to follow, we only 
consider the small-scale fading (i.e., Kb^k = 1 for all b and k). These setups can be considered to 
correspond to the worst-case scenario where all users are at the cell edge. A more realistic channel 
model where large-scale fading is taken into account is investigated in Fig. [5] for a massive MISO 
system. All noise variances are taken as unity and the transmit power is normalized with respect 
to the noise variance. For the sake of simplicity, but without compromising generality, we take 
Ofc = 1 for all k. A;, ^ = I^^, and pb^k = P for all 6, k. Unless otherwise mentioned, the error 
vectors are assumed to lie in a hypersphere of radius p. We evaluate the performance of the 
three approaches in terms of the worst-case SINR (i.e., the objective obtained at the end of the 
bisection procedure when solving dD and (|22]) ) and the probability of exceeding the worst-case 
SINR which is referred to as PE from now on. For the simulation setup considered in this 
paper, the SDP-based approach is numerically found to produce precoding matrices close to 
rank-1 matrices. 

Fig. [2] plots the average worst-case SINRs (over 200 realizations of the nominal channels hf, fc) 
versus the radius of the uncertainty sets, p, for all approaches. In this simulation setup, we only 
consider the small-scale fading (i.e., Kf, ^ = 1 for all b and k). The number of users in each BS 



is 2, i.e., there are A' = 4 users in total. It is seen in Fig. |2(a)| that the SOCP-based solution 
gives the worst-case SINR close to that of the SDP-based approach and slightly higher than that 
of [m when p' = p. As mentioned earlier, by taking p' as a design parameter, we can achieve 
a trade-off between the worst-case SINR and the resulting PE. It can be seen from Fig. [2] that 
the worst-case SINR of the SOCP-based solution is improved when we take p' = p/2.5, but this 
implies reduced PE as indicated in table |2(b)"j The values of PE given in table |2(b)] are obtained 
with 10^ sets of channel errors that are uniformly distributed in the ball of p using the toolbox 
of ETl . As expected, the non-robust approach delivers the maximum SINR and virtually zero 
PE in all cases. 
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(b) Variation of PE with p for the proposed SOCP-based robust 
design. 

Fig. 2. Average worst-case SINR versus p for different approaches where channel uncertainties are bounded by /2-nomi. The 
value of power for both BSs is taken as 5 dB. The number of transmit antennas at each BS is T = 8. The total number of users 
is _?sr = 4 (2 users per base station). 



Although the SDP formulations can offer better worst-case SENR, they are not practically 
useful for large-scale antenna systems especially from the complexity perspective. Another 
disadvantage of the SDP approach is its inability to handle various uncertainty sets and limited 
choice of solvers compared to those for SOCP based solutions. The flexibility in choosing a 
solver is important because a general purpose convex programing solver may not be efficient 
for all problems. We compare the simulation time of the SDP and SOCP-based robust designs 



TABLE I 

The average run time (in seconds) versus the number of transmit antennas, T, at each BS for the robust 

DESIGNS. The number of BSs is B = 2, EACH serving 10 users. The bisection procedure terminates when the 

difference between the objective values of two bisection steps, e < 10~^ . 
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SOCP-based design (GUROBI) (sec) 


1.08 
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90.72 



using YALMIP II28B with two widely used conic programming solvers (SeDuMi [291 and SDPT3 
Il30l ). Note that the proposed SOCP-based method allows us to make use of GUROBI |[3ll as 
a solver as well which is claimed to be very efficient for detecting feasibility of large-scale 
SOCPs. For the robust SOCP-based design, we use the simplified representation in (|23T |. In 
Table H we show the average run time (in seconds) of all robust approaches as a function 
of the number of transmit antennas, T, for solving the corresponding optimization problem. 
The bisection procedure terminates if feasibility is detected and the relative difference e of the 
objectives between two bisection steps is less than or equal to 10~^. The lower threshold of 
the bisection algorithm is set to 0, while the upper one is equal to the balanced SINR obtained 
from the non-robust design. The codes are executed on a 64-bit desktop that supports 8 Gbyte 
RAM and Intel CORE 17. For both solvers, it can be clearly seen that the SOCP-based design 
requires a lower run time and the difference is considerable as the number of transmit antennas 
T increases. This observation matches with the theory presented in the subsection on reduced 
complexity in Sec. IIII-B[ Moreover, we notice that the SDP-based methods are not capable of 
producing a solution when T > 50 due to lack of memory (denoted by a cross mark "x" in 
Table IB- When T > 300, SeDuMi and SDPT3 are not suitable solvers for the SOCP-based 
method since they are not able to produce a solution even after several hours (denoted by "— " 
in Table IJ)- We have observed that GUROBI is the most efficient solver for the SOCP-based 
method in particular for large-scale antenna array systems. 

In Fig. [3l the average worst-case balanced SINR (again over 200 realizations of the nominal 
channels h^^k) is plotted with the transmit power per BS, P, for different approaches. We note that 
the SOCP-based approach performs nearly as good as the SDP one. The reduced minimum SINR 
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Fig. 3. Comparison of the worst-case SINR of non-robust, SOCP-, SDP-based design, and the approach of (6] as a function 
of the transmit power of the BSs for p — 0.4. The performances of the SOCP-based method are shown for three cases where 
p' = p, p' = p/1.25 and p = p/2.5. 



of [[6l can be probably attributed to the fact that it is completely a conservative approximation 
of the robust counterpart. Recall that, in our proposed SOCP-based design, we can control the 
degree of conservatism of the design by finding proper value of p' . The values of PE for three 
approaches are also provided in Fig. |3] Further, being oblivious to channel error vectors, the non 
robust design delivers the best worst-case SINR. Nonetheless, as expected and seen previously, 
this comes at the cost of unacceptably low PE, i.e., PE ^ 0. It is found that the SOCP-based 
method gives PE = 0.99 for p' = p/1.25, while the approach of flU and the SDP-based solution 
both produce PE = 1.0. The value of PE for the SOCP-based design is reduced to 0.82 as we 
set p' = p/2.5. Interestingly, this decrease in PE is accompanied by a corresponding increase 
in the worst-case SINR, thereby providing a tradeoff between the two parameters. We note that 
the trend of values of PE is observed to be typical for the range of transmit power considered 
in Fig. [31 

In Fig. |4] we evaluate performance of robust beamforming for SINR balancing where the 
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(a) Variation of average worst-case SINR versus p for Zoo-norm. 
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(b) Variation of PE with p for Zoo -norm for the proposed SOCP- 
based robust approach. 

Fig. 4. Variation of average worst-case SINR versus p for Zoo -norm (i.e., "box" uncertainty). The value of power for both BSs 
has been taken as 5 dB.The number of transmit antennas at each BS is T = 8. The total number of users is K = 4. 



errors in elements of channel vectors are bounded within a (multi-dimensional) box of size p, 
i.e., |[vbfc]j| < p for all i. This is equivalent to saying that ||v(,fc||oo < P- For this case, we note 
that the SDP formulations and the approximations used in H are not applicable. The curves in 
Fig. m have been obtained by noting the fact that the dual of /oo-norm is /i-norm. It is seen that, 
with box uncertainty (Fig. lU, the worst-case SINR is lower than that in the case of ellipsoidal 
uncertainty for the same p. This can be explained as follows. We note that for a vector v, 
||v||i > ||v||2 > ||v||oo, which means that for the same p, the /oo-norm defines a smaller feasible 
set in (|22)) compared to the /2-norm. Thus the worst-case SINR for the /oo-norm uncertainty is 
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Fig. 5. Variation of average balanced SINR versus the number of transmit antennas per BS, T. For robust designs, channel 
uncertainties are are bounded in a ball of radius p. The value of the normalized power for both BSs is taken as 20 dB. The 
total number of users is K — 20 (10 users per cell). 



lower than that of the /2-nomi. However, when errors are uniformly distributed in a box, we 
note in the table given in Fig. |4(b)| a slight degradation in PE when p' = p. This stems from the 
fact that the proposed approach is not guaranteed to be safe as discussed earlier in the paper. 

Finally, in Fig. [5] we investigate how the balanced SINR of all users scales with the number 
of transmit antennas. In particular, we consider a system of two cells, each serving 10 users. The 
users are uniformly distributed in the cell and are not allowed to be closer to the BS by more 
than do = 100 meters [|32l . We also assume that the cell diameter (to a vertex) is 1000 meters. 
The large-scale fading coefficient is modeled as n^^k = /3fe,fc('^b,fc/'^o) '^ where (3b^k accounts for 
shadow fading assumed to follow a log-normal distribution with standard deviation (Jshadow^ ^ is 
the path loss exponent, and db^k is the distance between BS b and user k. In Fig. [5l we choose 
^shadow = 8 dB and i^ = 3.8 as in |[32ll . In Fig. [51 the performance of zero-forcing beamforming 
(ZF-BF) scheme is also included for comparison. In ZF-BF, multiuser interference for each user 
is forced to zero, i.e., hb.jini = for all j ^ i ll33l . Il34ll . In this way, problem dH) is simplified 



as 

maximize t (24) 

where 7" = {{t,uik)\hb^jmi = 0, Vz ^ j, \\Mb\\2 < y/Pb, V6, a-^ < ^{hb^^ki^k), V/c} which is 
convex. Consequently, problem (l24l) is jointly convex in t and m^. Using the null-space technique 
as devised in ll33l . [|34l . we can further remove the ZF constraints h^. jiiij = 0, Vi 7^ j in J^ 
without loss of optimality as follows. Let Mt = [hj hT ■ ■ ■ h7 hT • • ■ h7 1^ G 
(j-^(ft:-i)xT ^j^j^ Q^ ^ £Ty.{T-K+i) ^^ ^ matrix of orthogonal columns that span the null space of 

MfcO Then, to satisfy the ZF constraints, we can write m^ = Gfcihfc. The problem (|24] ) is thus 
further equivalent to 

maximize t (25) 

where hb^,fc = hb^^fcGfc and M;, = [mu,(i), ...,mu,(\Ut\)] and J" = {mfc|||(M5)||2 < y/Pi^yb, a-^ < 



otk 

3?(hfe^,,femfc), V/c}. The advantage of the ZF-BF scheme is that we can avoid the bisection 
procedure that must be carried out to solve (2) for optimal linear beamforming. However, even for 
perfect CSI, the performance of ZF-BF is still far away from that of optimal linear beamforming 
as shown in Fig. [51 For example, when T = 120, a gap of about 4 dB is observed between 
ZF-BF and general linear beamforming. These non-robust designs are sensitive to channel errors 
as shown by the fact that PE ?» for both these cases. The benefit of using large-scale antenna 
systems is seen in Fig. [51 when a gain of 13 dB is observed as the number of transmit antennas 
is increased from T = 8 to T = 120 for optimal linear beamforming with perfect CSI. A similar 
conclusion also applies to the robust designs when p is taken as 0.001, 0.005 and 0.01. We note 
that these values of p are comparable to the average channel gains of h^ ^ taking into account 
the effect of shadowing and path loss. Therefore, owing to the conservative nature of robust 
designs it is not possible to achieve nontrivial SINRs for higher values of p. 

V. Conclusion 

We have studied the design of beamformers that balance the SINR of users in a multicell 
downlink system in the presence of channel uncertainties. Norm bounded channel uncertainty 
model is used. As a first approach to solving the problem in this scenario, we present an S- 

^For the ZF-BF scheme to be feasible, i.e., dim(Gfe) > 0, we must have T > K. 



lemma based approximate solution in which the beamformers are obtained by solving an SDP 
in conjunction with bisection search. Later, by exploiting various properties of the functions 
involved in the problem, we present a solution in which robust beamformers are solutions to 
an SOCP-based formulation. We show that in addition to being capable of handling different 
uncertainty sets, the SOCP-based approximation exhibits a much reduced complexity solution. 
We have tested the performance of proposed approaches for the recently conceived massive 
antenna systems, and have determined that the SOCP approach outperforms the SDP-based 
solution from computational cost perspective. Finally, we have also shown that the reduced 
complexity SOCP-based approach yields a balanced SINR and the probability of achieving it 
which is comparable with the SDP approach. 
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Appendix 
Proof of Theorem [H 

We will follow the notation in (fT9l) and obtain a tractable version of this constraint by adapting 
the arguments developed in fl7]. Let Oi and O2 be the optimal solutions of 

max ci^^vi + a2 V2 (26a) 

subject to ||vi + V2II < "(9 (26b) 

vi>0,V2>0. (26c) 



and 



max ^max{[ai]j, [a2]j,0}[v3]j (27a) 

subject to llvsll < ?9. (27b) 



respectively. It is shown in IfTSl that Oi = O2. Now consider the following set of relations that 
can be obtained from (fT9] ). i.e., 



A(M,,z,,,,h,,,.)>-,„ min , ^ /2(M,,(5^,,)|[e,,fc].| + /2(M,,-<,)|[c|jJ,| 

\\\Sb,k\+\<i>b.M<p't,k ~^ I J 



(28a) 



H4'6,fclll<pU'^ I J 



(28b) 



lb,k f 

= „ "^^^ , Zl ^ "^^^ ( " ^2(M6, 5^), -/2(Mb, -5^^;,), 0) |[v36,3fc]i| 

l|v36,3fcll<Pb,fc -^ L 

(28c) 
where in (I28cl l we have employed the result in optimization problem formulation (1271) . and have 
slightly changed the representation in (fT9] ) by not specifying any particular norm in the constraint 



set. Now recalling the definition of the dual norm ?9||y ||* — max s^y of vector y, we obtain the 

||s||<'i9 

result stated in Theorem [1] for the constraint set of interest. 
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